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NOTATION 

particle  concentration  in  the  sample,  g  cm'"3 
particle  concentration  in  the  free  stream,  g  cm"3 
particle  diameter,  cm 

distance  from  the  inlet  to  the  outlet  cross  section  of  the 
collection  tube,  cm 

thickness  of  the  collection  tube  wall  at  the  outlet  cross 
section,  cm 

length  of  coaxial  boundary  tube,  cm 
radial  co-ordinate  of  particle  position,  cm 
radius  of  coaxial  boundary  tube,  cm 
radius  of  collection  tube  at  exit,  cm 
inlet  radius  of  cone,  cm 

radial  co-ordinate  of  particle  position  far  upstream,  cm 

far  upstream  radius  of  the  stream  tube  that  impinges  on  the 
collection  tube  circumference,  cm 

time,  seconds 

radial  component  of  local  fluid  velocity,  cm  sec-1 

axial  component  of  local  fluid  velocity,  cm  sec-1 

fluid  velocity  in  collection  tube,  cm  sec”1 

fluid  velocity  at  boundary  tube  entrance,  cm  sec”1 

fluid  velocity  at  collection  tube  exit,  cm  sec"1 

fluid  velocity  at  boundary  tube  exit,  cm  sec”1 

fluid  velocity  at  inlet  of  sampler,  cm  sec”1 

radial  component  of  local  particle  velocity,  cm  sec"1 

axial  component  of  local  particle  velocity,  cm  sec”1 

axial  co-ordinate  (origin  at  collection  tube  inlet)  of 
particle  position,  cm 


NOTATJON  (Cont'd) 


zQ  axial  co-ordinate  of  particle  far  upstream,  cm 

u  absolute  viscosity  of  fluid,  poise 

p  fluid  density,  g  cm'3 

o  particle  density,  g  cm*3 

1 1>  stream  function,  cm3  sec'1 


NOTATION  (Cont'd) 

spherical  particle  Reynolds  number  in  free  stream 

radial  component  of  local  fluid  velocity,  du/dr' 

axial  component  of  local  fluid  velocity,  du/d'z 

radial  component  of  local  particle  velocity,  dr/dt 

axial  component  of  local  particle  velocity,  dz/dt 

axial  co-ordinate  (origin  at  collection  tube  inlet)  of 
particle,  z/r 

axial  co-ordinate  of  particle  far  upstream,  zo/rc 

axial  co-ordinate  used  in  calculation  of  the  stream  function 
field,  z/rft 

ratio  of  collection  tube  radius  to  boundary  tube  radius,  r^/r^ 
length  of  coaxial  boundary  tube,  L/r^ 

distance  from  the  Inlet  to  the  outlet  cross  section  of  the 
collection  tube,  D/rA 

distance  from  inlet  of  boundary  to  inlet  of  collection  tube, 

3  -  Y 

time,  tllA/rc 

dimensionless  group  independent  of  particle  position,  ReQ2/K 


stream  function,  ip/^U^r^2 
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ABSTRACT 

Sampling  and  collection  efficiencies  are  calculated  for  a 
large-yolume  air  sampler  under  conditions  of  anlsoklnetlc  as  well  as 
Isokinetic  flow.  A  mathematical  model  developed  to  evaluate  a  tapered- 
tube  sampling  probe  was  modified  to  obtain  results  for  the  large- 
volume  sampler,  using  various  particle  sizes  and  flow  velocities.  These 
results  should  facilitate  the  prediction  or  correction  of  sampling  errors 
In  field  and  laboratory  experiments. 
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1 .  INTRODUCTION 

In  order  to  assess  the  effectiveness  of  a  specific  large- 
volume  air  sampler  (cyclone  scrubber),  the  Instrument's  ability  to  collect 
samples  of  finely-divided  particulate  matter  must  be  determined.  These 
samples  can  come  from  still  or  moving  alrstreams,  and  can  vary  both  In 
particle  size  distribution  and  In  concentration.  The  bio-sampler  under 
evaluation  consists  of  an  air  Inlet  cone  and  collection  unit,  and  Is  de¬ 
signed  to  operate  at  a  capacity  of  1000  litres  (air)  per  m*nute.  (It  Is 
described  fully  In  Suffleld  Technical  Note  No.  311). 

Sampling  from  streams  of  suspended  particulates  is  representa¬ 
tive  only  if  the  size  distribution  and  content  of  particles  In  the  sample 
are  Identical  to  those  of  particles  In  ambient  air  at  the  point  of  sampling. 
The  sampling  system  may  give  rise  to  three  different  types  of  error  (Vltols, 
1964)  due  to: 

(1)  parLlcles  falling  to  enter  the  sampling  cone  in  representa¬ 
tive  concentrations; 
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(2)  particles  being  deposited  between  the  air  inlet  cone  and 
the  collection  location;  and 

(3)  particles  being  shattered,  agregated  or  incompletely  re¬ 
tained  by  collection  devices. 

When  the  velocity  of  gas  entering  the  inlet  cone  is  exactly  the  same  as 
the  far-upstream  velocity  of  the  gas  ('isokinetic'  sampling),  particles 
will  enter  the  sampler  in  representative  concentrations.  Otherwise, 
errors  of  the  first  type  will  occur  as  the  result  of  anisokinetic  sampling. 

The  purpose  of  this  report  is  to  describe  the  modification  of  a 
mathematical  model  devised  formerly  for  calculating  the  error  due  to 
anisokineticlty  (Mellsen,  1979)  of  a  sampling  probe  developed  and  used  at 
DRES.  The  model,  previously  applied  to  a  straight,  tapered  tube  Is  herein 
adapted  to  the  funnel-shaped  inlet  cone  of  a  specific  large-volume  air 
sampler,  and  as  such,  calculates  the  sampling  and  collection  efficiencies 
produced  by  varying  upstream  gas  velocity  and  particle  size. 


2.  DEFINITION  OF  THE  PROBLEM 


As  explained  In  Suffleld  Technical  Paper  No.  499  (Mellsen,  1979), 
the  problem  of  finding  the  sampling  and  collection  efficiencies  is  one  of 
determining  the  values  of  the  upstream  particle  and  fluid  radii.  The  up¬ 
stream  particle  radius,  r  ,  is  defined  as  the  radius  of  the  limiting 

P  J»® 

particle  trajectory  envelope  which  encompasses  all  particles  (of  any  given 
diameter)  entering  the  sampler.  The  upstream  fluid  radius,  r,  ,  is  the 

J  I 

radius  of  the  stream  tube  impinging  on  the  outer  circumference  of  the  in¬ 
let  cone,  and  containing  the  total  volume  of  air  passing  through  the 
sampler.  The  sampling  efficiency,  proportional  to  the  areas  of  upstream 
particle  envelope  and  fluid  stream  tube,  can  then  be  calculated: 


C 

Co 


(Eq.  1) 


where  Co  is  the  upstream  particle  concentration  and  C  is  the  particle 
concentration  in  the  sample;  the  collection  efficiency  is  given  by: 


Em  = 


(Eq.  2) 
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where  rc  is  the  radius  at  the  inlet  of  the  cone. 

Inertial  and  drag  forces  may  cause  particles  flowing  far  up¬ 
stream  of  the  collection  inlet  to  deviate  from  stream  lines  on  arriving 
at  the  cone,  where  the  fluid  velocity  may  be  changing  markedly.  Thus, 
in  obtaining  the  true  free  stream  concentration  of  particles  and  the 
sampling  efficiency,  the  two  different  values  of  upstream  particle 
radius  and  upstream  fluid  radius  must  be  known.  When  the  free  stream 
velocity,  UA,  Is  less  than  the  sampler  inlet  velocity,  jl.e.  <  1  , 

some  particles  originally  inside  the  limiting  stream  tube  will  pass  out¬ 
side  the  sampler,  whereas  for  UA  , ,  some  particles  originally  outside 

u- 

the  stream  tube  will  be  drawn  into  the  sampler. 

3.  DESCRIPTION  OF  THE  SAMPLER 

The  part  of  the  large-volume  air  sampler  which  determines  stream 
function  values  and  hence,  affects  sampling  and  collection,  Is  the  air  In¬ 
let  cone  (Figure  1).  With  an  inlet  radius  of  2  1/2  inches,  the  cone  con¬ 
verges  to  a  straight  tube  of  Inside  radius  3/8  Inch,  through  a  funnel 
shaped  by  the  Intersection  of  three  circular  arcs.  The  entire  Inlet  cone 
Is  6  Inches  long,  the  converging  section  being  4  Inches  and  the  straight 
tube,  therefore,  2  Inches.  The  wall  of  the  cone  Is  1/16  Inch  thick,  but 
although  this  was  taken  Into  account  In  the  calculation  of  the  velocity 
Uc>  the  wall  thickness  was  neglected  In  the  computations  leading  to  the 
array  of  stream  function  values.  Since  a  grid  unit  in  the  array  repre¬ 
sents  1/8  inch,  the  cone  wall  thickness  of  1/16  Inch  would  have  little 
effect  on  stream  function  values,  but  would  make  computing  procedures  un¬ 
necessarily  complicated, 

The  three  circular  arcs  defining  the  shape  of  the  inlet  cone 
are  (numbers  in  Inches): 

(X!  -  2  5/16)^  +  (Y i  -  3  5/8P  «  (2«*)a  at  inlet  (Eq.  3) 

(X;,  -  4  3/8 )?  +  (Y?  -  6  16/32)^  •  (6)a  In  middle  (Eq.  4) 
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(X:i  -  4 )■'  +  (Y3  -  3)?  =  (2  1 9/32 )2  just  before  (Eq.  5) 

straight  tube 

Newton's  Method  was  used  to  determine  the  two  intersection  points 
(between  Equations  3  and  4,  and  Equations  4  and  5),  with  initial  values 
for  the  iterative  technique  found  by  inspection  of  a  drawing  of  the  curve. 


4.  EQUATIONS  OF  MOTION 

The  equations  of  motion  were  established  in  a  previous  report 
(Mellsen,  1979),  but  are  included  here  for  completeness. 


The  motion  of  an  Individual  particle  has  been  shown  (Vitols,  1964 
and  Batchelor,  1956)  to  be  determined  by  the  following  ordinary  differential 
equations : 


dVp.  CpRe('ty  -  v--) 
cfT  3  24K 


(Eq.  6) 


dv2  CpRe  ( ijT-  -  vz-) 
dT'  "  24K 

where  Re  =  ReQ I (u'r  -  v--)2  +  (u--  -  v-)2]'5 

o  d2U. 

K  =  yg,  r-  particle  inertia  parameter 
U.dn 

Re„  =  free  stream  Reynolds  number 

0  |l 


(Eq.  7) 
(Eq.  8) 

(Eq.  9) 

(Eq.  10) 


The  symbols  are  defined  in  the  notation  section  near  the  front 
of  this  report  and  the  basic  geometry  of  the  flow  system  is  illustrated 
in  Figure  2. 


Several  assumptions  are  inherent  in  the  use  of  Egs.  6  and  7 
for  calculating  the  collection  and  sampling  efficiencies  due  to  a  stream 
of  particles,  including: 

(a)  uniform  particle  distribution; 

(b)  no  gravitational  or  electrostatic  forces  of  consequence; 

(c)  monodisperse  spherical  particles  with  diameter  very  small 
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in  relation  to  the  inlet  diameter  of  the  sampler;  and 

(d)  free  stream  flow  that  is  steady,  imcompressible  and 
irrotational . 

The  drag  coefficient  is  a  function  of  Reynolds  number  and  is 
available  in  the  form  of  definitive  empirical  equations  (Davies,  ’145). 
These  equations  are  stated  as  follows: 

C  Re 2 

Re  =  JL —  -  2.3363  x  1 0~‘*  (CnRe2)2  +  2.0154  x  10"6  (CnRe2)3  - 
24  u  u 

6.9105  x  10'9  (CDRe2)4  (Eq.  11) 

for  Re  <  4  or  C^Re2  <  140 

logi0  Re  =  "  1.29536  +  9.86  x  10"1  (1og10CDRe2)  -  4.6677  x  10"2 

(log:  cfDRe2)2  +  1.1235  x  10‘3  (logi  <£DRe2)3  (Eq.  12) 

for  3  <  Re  <  10“  or  C^Re2  <  4.5  x  107 

5 .  AIR  FLOW  FIELD  EQUATIONS 

These  equations  were  stated  and  explained  in  an  earlier  report 
(Mel  1  sen,  1979),  but  are  again  shown  for  the  sake  of  thoroughness. 

The  equations  of  fluid  velocity  were  derived  from  the  stream 
function  for  ideal  flow  over  and  through  the  sampler.  To  solve  the  prob¬ 
lem,  an  outer  boundary  was  used  around  the  collection  cone  in  the  form  of 
a  coaxial  tube  of  radius  (Figure  3),  which  was  chosen  large  enough  so 
that  the  effect  of  the  boundary  tube  on  flow  in  the  proximity  of  the 
sampler  is  negligible.  The  collection  cone  was  inserted  a  distance  D 
into  the  downstream  end  of  the  boundary  tube.  Since  the  flow  is  axi- 
symmetric  only  a  radial  plane  containing  Doth  tubes  had  to  be  considered. 

The  fluid  enters  the  boundary  tube  with  steady  velocity  Uft, 
and  separates  into  a  central  stream  with  velocity  at  the  exit  and 
at  the  entrance  of  the  sampler,  and  an  annular  stream,  with  velocity  U^, 
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at  the  downstream  end  of  the  boundary  tube.  The  axial  velocities  U^, 

UD,  Ur  and  U.  are  uniform.  Also,  there  is  no  radial  flow  at  the  end 

D  L  1 

cross  sections. 

The  boundary  conditions  on  the  flow  can  now  be  completely 
specified  so  that  the  flow  field  can  be  obtained  by  solution  of  the 
equation  of  the  stream  function. 

The  axially  symmetric  stream  function  i^(r%a)  (Batchelor,  1967) 

satisfies: 


-  \  +  -l-t  =0  (Eq.  13) 

3 r?  r  3r  3z2 

The  two  velocity  components  (Figure  2)  are  given  by: 


u  =  -1- 
z  r  3r 


u  =  -  1 
r  r  32 


(Eq.  14) 
(Eq.  15) 

When  UA  and  Ug  are  specified,  continuity  gives  Uc  as  follows: 

(Eq.  16) 


IL  =  UA  fl2UB 

c  —  +nrrT1 


1  - 
r 


rB  P 


where  «  = 


B 


(Eq.  17) 


and  h  is  the  thickness  of  the  collection  tube  wall 


For  uniform  velocity  profiles,  the  stream  function  is  of  the 


form: 


ip  =  H-ur2 


(Eq.  18) 


To  allow  for  greater  generality,  the  stream  function  and  the 
geometric  variables  are  restated  in  the  following  dimensionless  form: 
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=  t 

' 

(Eq.  19) 

II 

(Eq.  20) 

.  z 
"  rA 

(Eq.  21) 

.  L 
"  rA 

(Eq.  22) 

_  D 

(Eq.  23) 

"  rA 

J.C 

II 

(Eq.  24) 

A 


The  boundary  values  for  the  stream  function  and  the  geometric 
configuration  in  terms  of  the  dimensionless  variables  are  shown  in  Figure 


The  axially  symmetric  stream  function  equation  (Figure  13) 


becomes: 


I2^  1  + 

3R2  ‘  R  »R  3Z2 


=  0 


(Eq.  25) 


6 .  DISCRETIZATION  SCHEME  FOR  THE  AIR  FLOW  FIELD 

The  equation  for  the  axially  symmetric  stream  function  (Eq.  25) 
is  discretized  as  follows: 


Vi.j  -  2vij  +  Vi.j 

ar2 


Vi.J  ~  Vi  J  + 

21  aR2 


!±±LL?!±±  -  0  (Eq.  26) 

AZ2 

where  1  and  j  are  the  grid  point  numbers  In  the  R  and  Z  directions 
respectively  (Figure  5).  Eq.  26  can  be  rearranged  to  give  a  simple 
equation  by  choosing  a  square  grid  so  that  AR  and  aZ  are  equal.  The 
resulting  equation,  which  is  suitable  for  Gauss-Seldel  Iteration 
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(Carnahan  et  al ,  1  969),  is  given  as  follows: 


Vl.J  *  Vl.J  *  T1,J-1 


'  1  .j+1 


(Eq.  27) 


Eq.  27  can  be  applied  to  all  Interior  points,  which  are  defined  as  points 
for  which  the  nearest  boundary  is  at  least  one  grid  unit  away. 


In  dealing  with  points  on  or  surrounding  the  boundary  described 
by  the  sampling  cone  (for  which  the  nearest  boundary  in  either  the  horizon 
tal  or  vertical  direction  is  less  than  one  grid  square  away),  a  Taylor 
series  expansion  was  used  (Carnahan  et  al,  1969)  and  the  following  finite 
difference  equations  derived.  (The  first  applies  to  points  below  or  to 
the  left  of  the  curved  boundary,  and  the  second,  to  points  above  or  to  the 
right  of  the  boundary.) 


ab 

a+b 


_  ab 
1 ,  j  "  a+b 


% 

+  a*(a+T7  +  bflj+TT  " 
.  yH 

+  a(a+l )  +  b(b+l T  " 


(Eq.  28) 
(Eq.  29) 


where  a  is  the  vertical  distance  (0  <  a  i  1 )  to  ^  and  4*y  represents  (for 
points  below  the  curve)  either  the  boundary  Y-value  (If  the  boundary  lies 
between  ^  j  and  H,i+1  j)  or  the  adjacent  f-value  j).  (For  points 

above  the  curve,  f  takes  either  the  boundary  value  or  the  value  of  T  ,  . 

V  » ""  ‘  >  J 

In  the  horizontal  direction,  b  is  similarly  defined  as  the  distance 

(0  c  b  L  1)  to  iu,  and  yu  is  the  closer  of  the  two  i-values,  the  boundary 

n  n  „  „ 

value  and  the  adjacent  value  (  i , j+1  for  points  below  the  curve  and  1,j-l 
for  points  above  it) . 


The  grid  size  was  chosen  from  past  experience  (Mellsen,  1979)  so 
that  each  grid  unit  (both  horizontally  and  vertically)  represents  1/8  Inch. 
Transferred  to  the  grid  (Figure  5),  the  straight  tube  radius,  r^,  then 
corresponds  to  3  units,  the  inlet  radius,  rc>  is  20  units,  the  boundary 
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tube  radius,  r^,  is  120  units  (to  be  located  a  distance  of  five  times 
the  inlet  radius  outward  from  the  edge  of  the  cone),  the  complete  length 
of  the  inlet  cone,  y,  is  48  units,  and  the  distance  to  the  upstream  end 
of  the  boundary  tube,  A,  is  152  units  (so  as  to  be  more  than  seven  inlet 
radii  upstream  of  the  collection  inlet).  Specifying  the  boundary  tube 
radius  and  the  distance  to  the  upstream  boundary  in  this  way  ensures  that 
the  behaviour  of  the  flow  be  as  if  the  inlet  cone  were  situated  in  free 
space  and  the  particles  coming  from  such  a  distance  upstream  as  not  to  be 
affected  by  the  cone. 

The  stream  function  was  obtained  by  Gauss-Seldel  iteration  using 
Equation  27,  28  and  29.  The  boundary  conditions  were  set  initially  at  the 
centerline,  at  the  boundary  tube  wall  and  inlet,  and  at  the  outlet,  and 
held  constant  throughout  the  iterative  procedure.  Any  point  not  falling 
on  either  one  of  these  boundaries  or  the  wall  of  the  inlet  cone  was 
initialized  to  zero.  A  Fortran  program  (listed  In  Appendix  A)  was  written 
to  perform  the  calculations  on  an  IBM  370  computer. 

A  special  routine  (adapted  from  Carnahan  et  al,  1969)  to  handle 
points  near  the  curved  wall  of  the  inlet  cone  had  to  be  Incorporated  Into 
the  Fortran  program.  This  routine  first  labels  points  as  being  one  of  four 
types  (see  Figure  6)  by  finding  the  highest  point,  UMAX  (the  maximum  within 
the  boundary),  for  each  row,  I,  and  classifying  points  according  to  the 
horizontal  and  vertical  distances  to  the  curve  (B  and  A,  respectively): 

B  =  AK  -VRK~-~f(I-lT  -  BK~p-'  -  (J-l  )  (Eq.  30) 

A  =  BK  ("(J-l  j~AKp;  -  (1-1)  (Eq.  31) 

where  I  and  J  are  the  coordinates  of  the  point,  and  AK,  BK  and  RK  assume 
the  values  of  the  a,  a  and  r  In  the  equation, 

(X  -  a)2  +  (Y  -  p)2  =  r2  (Eq.  32) 

from  the  particular  circular  arc  (Eqs.  3,  4  and  5)  defining  the  curve 
at  that  point.  The  distance,  A,  is  then  found  for  every  point  in  each 
row,  starting  at  JMAX  and  decreasing  along  the  row  until  an  Interior 
point  is  reached,  and  the  procedure  Is  repeated  on  the  right  side  of 
the  curve,  using  JMIN(I)  (the  minimum  above  the  boundary,  neglecting 
the  wall  thickness  of  the  cone): 
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JMIN(I)  =  JMAX(I)  +  1  (Eq.  33) 

and  continuing  until  the  upper  interior  point  is  reached.  The  horizontal 
and  vertical  distances  are  now  defined  by; 

BO  =  1  -  B  (Eq.  34) 

AO  =  (1-1)  -  [BK  -  ''v/RK?  -“(U-AK)?]  (Eq.  35) 


The  coefficients  of  the  v-values  in  Equations  28  and  29  are  then  cal¬ 
culated  using  A  and  B  (for  Eq.  28)  or  AO  and  BO  (as  a  and  b  in  Eq.  29). 

Type  IV  points  are  assigned  the  boundary  value  and  held  fixed 
through  the  program.  For  the  other  types  (I,  II  and  III),  the  values  of 
and  can  then  be  determined  and  the  iteration  performed  according  to 
Eq.  28  (for  points  below  the  curve)  or  Eq.  29  (for  points  above  the  curve) 
For  example,  for  a  Type  II  point,  j,  below  the  curve,  Tv  would  assume 
the  value  of  the  boundary  and  v^,  the  value  of  Vj*  while  If  j  were 

above  the  curve,  vv  would  again  assume  the  boundary  value,  but  would 

become  vi  .  j . 

7 •  SOLUTION  OF  THE  EQUATIONS  OF  MOTION 

From  Section  2  of  this  report,  as  In  a  prior  paper  (Mellsen, 
1979),  the  problem  is  to  find  the  upstream  particle  and  fluid  radii, 
rp  m  and  r$  respectively,  In  order  to  calculate  the  sampling  and 
collection  efficiencies.  In  the  same  dimensionless  form  of  Equations 
6  and  7,  the  value  of  rp  ^  (notation)  was  found  by  an  Iterative  procedure 
called  the  half  interval  method  (Carnahan  et  al ,  1969).  The  value  of 
rp  a>  for  a  critical  particle  was  estimated  far  upstream,  the  path 
followed  to  the  plane  of  the  cone  inlet  and  the  miss  distance  (from  the 
edge  of  the  inlet)  calculated.  Next,  the  aforementioned  half  interval 
method  was  applied  to  determine  a  better  Initial  estimate,  the  path  again 
followed  to  the  plane  of  the  inlet,  and  another  miss  distance  calculated. 
This  was  repeated  several  times  until  sufficient  accuracy  was  achieved. 

The  Initial  upstream  position  in  a  plane  perpendicular  to  the  flow 
direction  was  located  far  enough  from  the  Inlet  so  that  free  stream 
conditions  would  prevail.  A  distance  of  seven  inlet  radii  upstream  of 
the  inlet  was  considered  adequate  on  the  basis  of  the  five  Inlet  radii 
serving  the  case  of  strai  ,.it  tube  sampling  (Batchelor,  1956). 
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The  path  of  an  individual  particle  was  determined  step-by-step 

by  applying  a  fourth  order  Runge-Kutta  method  (Carnahan  et  al ,  1969)  to 

the  equations  of  motion  (Eqs.  6  and  7).  The  values  of  Re  and  K  in  these 

equations  were  easily  found  for  each  new  step  by  direct  substitution  of 

previously  determined  values  Into  Eqs.  8,  9  and  10,  but  the  value  of  CpRe 

in  Eqs.  6  and  7  had  to  be  calculated  in  each  step  by  numerical  solution 

of  the  definitive  empirical  equations  (Eqs.  11  and  12).  This  was  done 

using  Newton's  method  (Carnahan  et  al,  1969)  for  finding  the  zero  of  a 

function.  The  values  of  u-  and  uv  were  calculated  In  each  step  from  the 

r  z 

stream  function  field  as  follows: 


r  4(1-1 ) (AR)2 

u-  .  litlbLl hzLd 

2  4(1-1 ) (AR)2 


(Eq.  36) 
(Eq.  37) 


where  1  and  j  define  the  grid  point  of  the  particle  position,  Since  the 
inlet  radius  of  the  sampler  was  chosen  to  be  20  grid  units,  these  are 
given  by: 


1=1+  20r  (Eq.  38) 

J  »  JQ  +  20 (Z  -  ZQ)  (Eq.  39) 


where  j  and  zQ  are  the  starting  point  values  of  j  and  z.  The  values  of 

1  and  j  obtained  from  Eqs.  38  and  39  were  rounded  off  to  the  nearest 

lower  integer  value  in  each  calculation.  The  value  of  r„  was  obtained 

s 

directly  from  the  stream  function  by: 


(Eq.  40) 


calculated  at  the  lowest  value  of  1  satisfying: 


^.j0  >  s'ic.Jc  (Eq-  41 ) 

where  ic  and  jc  define  the  grid  point  at  the  edge  of  the  collection  cone 
inlet.  The  calculations  to  obtain  the  solutions  were  done  with  an 
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IBM  370  Computer  by  means  of  a  Fortran  program,  the  listing  of  which  is 

shown  in  Appendix  R.  The  sampling  and  collection  efficiencies  given  by 

Eqs .  1  and  2  were  also  obtained  by  this  program  after  the  values  of 

r  and  r  had  been  calculated. 

P.'”  s,- 

8.  RESULTS 

Method  of  Analysis 

A  stream  function  array  was  computed  for  each  of  the  following 

ratios  of  ^B.  400,  400,  400,  400  and  400.  Because  of  the  funnel  shape 

UA-  T  3  '9  Pf 


of  the  sampler,  tapering  from  an  inlet  radius  of  2%  Inches  to  a  straight- 
tube  radius  of  3/8  inch,  a  velocity  of  Ug  =  400  implies  an  Inlet  velocity 

M  A  +  4  SW*  l\ 

A 1 


of  9. 


This  means  that  the  sampling  velocity  ratios  (U 


are  1/9,  1/3,  1 , 


3  and  6.  When  the  sampler  operates  at  its  design  capacity  of  1000  Ji/min, 
the  values  of  Ug  and  then  become  Ug  =  5847.482  cm/s  and  =  131.5683 


cm/s,  so  that  UA  varies  from  14.62  cm/s /U 


n 

K 


\ 


to  789.4  cm/s 


A  broad  range  of  particle  sizes,  of  diameters,  6,  10,  20,  50,  100, 
200  and  500  microns,  composing  monodisperse  fields,  was  analyzed  for  each 
stream  function  array.  Results  were  tabulated  and  plotted  in  graphs  of 

r 

sampling  efficiency  versus  inertia  parameter  (Figure  7;  g-  vs  log  K), 

collection  efficiency  versus  Inertia  parameter  (Figure  8;  Urn  vs  log  K),  and 

sampling  efficiency  versus  sampling  velocity  ratio  (Figure  9;  C  vs  ^A). 

Co 


Discussion 

The  validity  of  the  mathematical  model  has  been  discussed  and 
reported  (Mellsen,  1979).  Results  of  the  present  work  (Table  of  Results) 
show  the  sampling  efficiency  of  a  uniform  field  of  20  micron  spherical 
particles.  In  a  wind  that  Is  six  times  the  sampling  inlet  velocity,  to  be 
in  error  by  over  30%  =  1.306).  Since  smaller  particles  are  carried 

more  readily  with  the  air  stream,  the  sampling  of  small  particles  gives 
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rise  to  smaller  error.  The  exact  errors  for  very  small  particles  cannot 
be  determined  by  this  model  because  computing  errors  increase  with  de¬ 
creasing  particle  size  (Mellsen,  1979),  the  reason  being  that  as  particle 
size  decreases,  a  larger  numbe1'  of  calculations  is  required. 


The  case  of  Isokinetic  sampling,  where  free  stream  velocity 


matches  inlet  velocity 


[UA 


- 1 


should  be  characterized  by  both 


sampling  and  collection  efficiencies  equal  to  1 


C 

Co 


1 .  Em  =  1 


This 


is  displayed  quite  well  by  the  predicted  values  of  the  model.  For  example 
(Table  of  Results),  a  6  micron-particle  field  Indicates  an  error  of  only 


1.7% 


C 

Co 


1.017  ,  and  a  100  micron-particle  field,  of  2.5% 


c.. 

Co 


.9748 


Although  the  model  cannot  be  used  for  the  prediction  of 
efficiencies  in  completely  still  air,  very  low  free  stream  velocities  can 
be  handled.  The  lowest  free  stream  velocity  currently  tested  and  plotted 
is  14.62  cm/s,  but  if  desired,  lower  velocities  might  be  tried.  The 
effect  of  varying  free  stream  velocity  while  keeping  the  sampling  velocity 
constant  is  clearly  Illustrated  in  Figure  9  for  selected  particle  sizes. 


9.  CONCLUSIONS 

The  effect  of  anlsokinetlclty  on  sampling  with  the  DRES  large- 
volume  air  sampler  Is  sufficient  to  produce  significant  errors  in  sampling 
and  collection  efficiencies.  A  mathematical  model,  formerly  applied  to  a 
straight,  tapered-tube  sampling  probe,  was  modified  to  be  applied  to  the 
specific  large-volume  sampler  developed  at  DRES.  The  results  from  the 
model  can  serve  to  predict  the  magnitude  of  sampling  errors.  Several 
free  stream  velocities  for  a  fixed  sampling  rate  were  evaluated  with  a 
number  of  monodisperse  fields  of  suspended  particles.  Therefore,  if  the 
wind  velocity  and  particle  size  and  density  are  known,  the  results  of  this 
model  can  be  used  in  correcting  measured  samples. 
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NOTES 

1 .  CONTOUR  DIMENSIONS  ARE  APPROXIMATE.  TRANSITION 
AND  POINTS  OF  TANOENCY  SHOULD  BE  SMOOTH. 


FIGURE  1 :  DESIGN  DRAWING  OF  LARGE  VOLUME  AIR  SAMPLER 
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FIGURE  5:  DIMENSIONS  OF  DISCRETIZATION  GRID  FOR  AIR  FLOW  FIELD 
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FIGURE  8:  TYPES  OF  BOUNDARY  POINTS 


UNCLASSIFIED 


UNCLASSIFIED 


21 


k 


i 


°3  N0I1VM1N33N03  MV1U18  33Ud 

0  NOIlVUiNaSNOS  31dWV8 


UNCLASSIFIED 


i 


\ 


FIGURE  7:  EFFECT  OF  VELOCITY  RATIO  ON  SAMPLING  EFFICIENCY 


FIGURE  8:  EFFECT  OF  VELOCITY  RATIO  ON  COLLECTION  EFFICIENCY 
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TABLE  I  RESULTS 


Ub 

Pa 

•Ja 

UA  (crn/s] 

1  d  (cm) 

K 

C 

co 

Em 

400 

1 

T 

9 

14.62 

.05 

1.787 

.3754 

2.749 

.02 

.2860 

.6708 

4.912 

.01 

.07149 

.8487 

6.215 

.005 

.01787 

.9358 

6.853 

.002 

.002860 

.9550 

6.993 

.001 

.0007149 

.9592 

7.024 

.0006 

-- 

-- 

-- 

400 

1 

1“ 

7 

43.86 

.05 

5.362 

.4645 

1.243 

.02 

.8579 

.6418 

1.717 

.01 

.2145 

.8014 

2.144 

.005 

.05362 

.9223 

2.486 

.002 

.008579 

.9681 

2.590 

.001 

.002145 

.9773 

2.615 

.0006 

.0007721 

.9804 

2.623 

400 

~r 

1 

131.6 

.05 

16.09 

.8950 

1.004 

.02 

2.574 

.9286 

1.042 

.01 

.6434 

.9748 

1.094 

.005 

.1609 

1.007 

1.130 

.002 

.02574 

1.015 

1.139 

.001 

.006434 

1.016 

1.140 

.0006 

.002316 

1.017 

1.141 

400 

~J7 

3 

394.7 

.05 

48.26 

1.628 

.9822 

.02 

7.721 

1.595 

.9624 

.01 

1.930 

1.516 

.9142 

.005 

.4826 

1.356 

.8182 

.002 

.07721 

1.161 

.7005 

.001 

.01930 

1.072 

.6466 

.0006 

.006949 

1.063 

.6413 

400 

■"5T 

6 

789.4 

.05 

96.51 

2.071 

.9805 

.02 

15.44 

2.029 

.9608 

.01 

3.861 

1.930 

.9137 

.005 

.9651 

1.699 

.8045 

.002 

.1544 

1.306 

.6185 

.001 

.03861 

1.154 

.5464 

.0006 

.01390 

1.102 

.5217 

UB  " 

5847.482 

cm/s, 

■  131.5683  cm/s 
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APPENDIX  A 

COMPUTER  PROGRAM  FOR  CALCULATING  THE  STREAM  FUNCTION 
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#Uir*T,  W  /RF/WFW8 ,  NP7I .  US,  WKiWZ.  NPR 
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i, 32), 0(20, 32), 

IS^0RAj;ui^6^NMS;NPZB.UB.RA,NPZ.NPR 
R£ AO  AND  CHECK  INPUT  PARAMETERS 

!it?l  !i?if  li  St  .*  s?'?*grsL  ,^7^^ 

calculate  and  write  0IMEN8I0NLESS  parameters 

ALPfU»RB/RA 


AKBAMWI"  - - “ 

2*UB ) / (!,»(( (RStO . 5 ) /RB) » ALPHA )**2 ) 

k«u^/yAJ*  ^ 
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NPZ*N2tl 


NP^SnIJi 


J.OAT  (NR) 


NPZBBNZB+l 
XRB*ALPMA*FL0AT(NR) 

NRB*IF I X ( XRB  ♦  0.1) 

NPRB*NRBAl 
NPRB1 »NPRBf 1 
NPRC*NRC+1 

Vo  IN  F1LE 

1.1,7 


IF ( INDEX  ) 
1  CONTINUE 


ESTABLISH  INITIAL  GUESSES  FOR  STREAM  FUNCTION 
AND  BET  BOUNDARY  CONDITIONS  ON  CENTRE  LINE 
AND  INLET  OF  OUTSIDE  FIRE 


I.ET 
■  1  .NR 


Rl*0£LR)**a! 


DO  2  1*1. N 

MifiJJUKo 

SET  BOUNDARY  CONDITION  AT  OUTLET  OF  INBXDE  PIPE 
PSI  ( I .  fiPZ )  ■  URAT*(RI«OELR)«*2 

SET  BOUNDARY  CONDITION  AT  OUTLET  OF  OUTSIDE  PIPE 
DO  «  I  *  NPRBl.NR 

mmuiiu  •UCRA*(1.0*(RI*DELR)**2) 

SET  BOUNDARY  CONDITION  AT  NECK  OF  INSIDE  PIPE 

NilHlHtyWM 

set  Boundary  condition  at  wall  of  outside  pipe 

DO  6  J*1  *  NR!! 

psknpr.jM.o 

set  boundary  condition  at  funnel  wall,1 


I 

7  H«N 


N*NpS-Np| 


ga, 

6-2 


HKiyi.RT  AI.8I.R1.A2.B2.R2.A3.B3.RS 
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nm 
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Sr  |i  ERs 
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FORMATS  FOR  OUTPUT  STATEMENTS 
200  FORMA T(6SH1 STREAM  FUNCTION  FOR  FLOP  IN  TWO  CO 
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32 


35 
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SUBROUTINE  F OH  DETERMINING  BOUNDARY  POINTS  AND  INTERCEPTS 
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MP 1  «M  + 1 
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JM  IN ( 1 ) *Nt2 
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14 

15 
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RKNR2 
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■  AO*BO/ ( AQt BO ) 
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locate  boundary  points  op  type  2 

DO  II  lal.H 
FIMHI-I 

jjp'/jHAX(l"l)»J)3<if  10*10 

innshil  il:U:lt 

AK  ■  A  t 


19 

20 
21 


RMRl 
SOTO  21 
AR»A2r 

«kbh| 

R*bR2 

goto  ai 

AKaAS 
BKBbJ 
HRaRl 


ABdK-SORr (HHBB2-(PJH1«AK)**2)«PIM1 
lTYPE(IiJ)aa 


1* 

15 

iH 

25 

2d 


JaJMAX (MP1 )•! 

ai 

If  iffs|:?HI  ii:ii:il 


i-n-j)  ro.ro. if 


-(pjHt-AK)**an 


rj,:}|s!:^o9(i6ii.) 

h!: !j)Si ,/?,/( ao+1 . i 

Ja,  tl 

ib3a 


^Ngl)ajHIN<MPl)  +  l 
C 

END 
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SUBROUTINE  ARCEPTIX. V , H 1 . P2 , ll , 12 . A 1 , AS) 

L’s|;Us!^II?*S5i:i:TiiiiT!i»STli:is;i5iTT«li:sii 
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i  yUa.SuSupiAAaSiiBAlSMa) 

a  c^nVinue 
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APPENDIX  B 


COMPUTER  PROGRAM  FOR  SOLVING  THE  EQUATION  OF  MOTION 


non  ooo  non  nnno  noon  nnnnnni 


/y^oTnm  jor-nrrwvi 

M  j m  fe- 

"U\ «y  iff!l:S 


'  »H»:K 


■(0.29) 

?2 ' OPT IHIZE (2 ) • 


NRlfI(S«5oIjG#LFT*GARITi  SIGNL.DT  AU»  NIBP.  NSBP*  NX 
^TAB^SH^ PHYSICAL  PROPERTIES  POR  CALCULATING  COLLECTION 

8$  II  HkfBB  RffiBITIB:  SB 

REAt)(5.*)DC,tp,RHa.mniA.iniuru».Tsr  . 


iMKPSI.rsmiMJieWSMTttiMJr 

.  MP-w  . 

S8kt.H"!8».?«?I,>8.*K4Io,0IMI»?iItMnT..» 


j.JI  INTVL 
fljsll)  INTVL 

IS  FREE  STRIAE  VELOCITY.  CN/BIC 


UA  «  UB/tIRAT 

3k»S?SSa*OP»*2«U^/(9.»XHUaDC) 


nitHe  Ti.loh  re r.  x*  rr.  dcv  up  .R  ho  ,  n  gh  a ,* n u/u  a . ur, hr 

ESTABLISH  GRID  STEP  SIZE 


NR  •  NPR  •  I 


HALF  INTERVAL  ITERATION  FOR  INITIAL  GA  VALUE 
09  21  ITERa|«  NX 

SET  AND  PNINT  INITIAL  CONDITIONS 


.  IP®' 

1  W  |eil-NA)s,7.a 

?  continue 


-G(l) )»A21**0.5 


dum 


TER.GRLFT.ORZER.ORRIT.TAU.BCI  J  #«(!>,. 0(1). Om.UZ. UR# 


CALL  ON- BUNGE  KUTTA  SUBROUTINE 


.ii 1  .li  2' .?  r 


'inn  no  o  non  non  non  non  non  non 


e  CONTI NUE 
H«Ht  1 

CALL  SUM22 ( 9 , G .0G  »  T  AU,  OT  AU ,  I  RUNG  .  M  ) 

IF  C  IRUMG-1 ) 10,9. 10 

9  REaREZ* ( (UR-G(2) )*«2*(UZ-G( l ))**?)* *0.5 
XCDREaCDHE (HE ) 

DG (  1 ) ■ ( ( XCDRE )/(29.0*XK))*(UZsG(l  )) 

DG(2)a( < XCDRE )/ <24. 0*XK )  )*(UR-G(2)) 

0G(3)»G( 1 ) 

DG(«)*G(2) 

GO  TO  8 

10  CONTINUE 

0 

CALCULATE  FLUIO  VELOCITY  AT  PARTICLE  POSITION 

I  •  1  ♦  1 F  IX(RC*G(4)  ) 

J  aJ0*lF IX(RC*(G(3)«G3ZtR) ) 

HUFLOATCI-1  ) 

UZMPSI  ( 1*1.  J)-PS1  ( 1-1 .  J))/(FDR8Q*RI) 

ur*  (ps  i  (i.  j-n-psi  n»  j*i  n/(F0RSQ*Ri) 

PRINT  SOLUTIONS 

IS  a  ITER/NIBP*N18P 
IF(I8-1TER)U.I3,11 

11  CONTINUE 
IF(ITER»1)12,13,12 

12  CONTINUE 

IF  1 1 TER-NX  >17# 1 i, 17 

13  CONTINUE 
NSTtP*NST£P* 1 
IF(NSTEP-N3BP)»7, 14, 17 

19  CONTINUE 
NSTEP-0 

TAR  ■  TAU  ♦  0.0001 

WRITE (6.209) TAW, G(1).G(2).G(3).G(4),UZ» UR. XCDRE 
INTEGRATE  ACROSS  ANOTHER  STEP  IF  REQUIRED 

17  HIT S*G( 3 ) 

IF (HITS)fl, 18, 16 

18  CONTINUE 

FIND  INTERVAL  HALF  WITH  THE  SIGN  CHANGE 

IF ( <G(9  1-1 .0)*SIGNL-0. 0)19, 19,20 

19  G9R 1 T *G9 ZE R 
GO  TO  21 

20  G9LF  T *G9ZER 

21  CONTINUE 

PRINT  SOLUTIONS  FOR  FINAL  VALUE  OF  TAU 
TAW  a  TAU  ♦  0.000) 

WHITE  U.20  9) TAW,G( 1 ) , G ( 2 ) , G ( 3 ) . G ( 9 ) . UZ . UR . XCDRE 

CALCULATE  THE  COLLECTION  EFFICIENCY 

WRI TE (6, 209 )  G9ZER 
EH  a  G4ZER**2 
WRITE(6,206)EH 

CALCUIATE  THE  3AHPLING  EFFICIENCY 

RSlNFaSGRT (PS  I  ( NPRC . NPZB ) ) *R A /RC 
PSIHT  ■  PSI  (NPRC, NPZB) 

DO  25  Ial.NPR 

IF (PSI (I, JO) -PSIHT) 25. 25. 29 

29  RSINF  a  F L 0 A T ( I -2) *SQR T ( PS IHT /PS I ( !- 1 , JO ) ) *DE LR*RA/RC 
GO  TO  28 

25  CONTINUE 
28  CONTINUE 

WRITE(6,210)  RSINF 
CHa(G9ZER/RSINF)**2 
WRITC(6.207J  CR 
REA0(5. *)N3TUP 
IF (NSTOP ) 1 . 30, 30 

30  STOP 
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FORMATS  FOR  OUTPUT  ST ATtMENTS 

200  FORMAT (  1  H  1  ,  J7X,  90HC0LLEC T ION  EFF 1 C I fc  N  C  Y  OF  A  CIRCULAR  TUBE/ 

1  1  HO  ) 

201  FORMA  T (  1 0HQU9LEF  a  ,P10.b/  10H  G9RIT  *  .FIO.b/lOH  SIGNL  a  , 

1  F  9 . 0/  10H  DTAU  =  .F10.6/  10H  NIBP  ■  .10/  10H  NSBP  *  .10/ 

2  10H  NX  s  ,101 

202  FORMAT (  10HOREZ  ■  .F12.7/10H  XK  s  ,E12.b/ 

1  10H  P  a  , E 1 0 . 0/ 

2  10H0DC  ■  ,  F10.5/10H  DP  a  .F10.7/10H  RHO  ■  .F10.6/ 

3  10H  SIGMA  a  ,  F10.6/10H  XMU  a  .F10.7/10H  UA  a  ,F10.0/ 

«  10H  UB  a  .F10.0/10H  ISR  »  ,15) 

203  FORMAT!  1 OHOI TER  a  ,13/  ,1  OH  GOLEF  a  ,F10.b/  10H  GOZER  *  , 

1  F10.6/  10M  GORIT  a  .FllO.b/  7H0  TAU.  11X.  OHG(l).  12X. 

2  OHG ( 2 ) ,  12X,  0HGC3).  12X,  OHG(O),  13X.  2MUZ,  10X.  2HUR  . 

3  1 2 X .  OHCOHE  / 

0  1  HO ,  F  7 . 0 ,  OF  lb. 6,  SFlb.O  ) 

200  FORMAT  (  1H  ,  F7.0,  OFlb.b,  3F16.0  ) 

205  F  OHM  A  T  C  ObHim  MOTION  OF  A  CRITICAL  PARTICLE  IS  GIVEN  BY  1 

206  FORMAT!  30H0THE  COLLECTION  EFFICIENCY  18  .  E10.0) 

207  FORMAT!  29H0THE  SAMPLING  EFFICIENCY  IS  ,£10.01 
loo  FORMAT!  28H1THE  PHYSICAL  PARAMETERS  ARE  1 


cun  r  unnp  i  i  in  9  r  r  *  n  #  nr  1 ui v 1  jr  1 u ■ m  4 

205  F  OHM A  T  !  ObH 1 T  HE  MOTION  OF  A  CRITICAL  PARTICLE  IS  GIVEN  BY  1 

206  FORMAT!  JOHOTHE  COLLECTION  EFFICIENCY  18  ,  E10.0) 

207  FORMAT!  29H0  THE  SAMPLING  EFFICIENCY  IS  ,£10.0) 

20fi  FORMAT!  28H1THE  PHYSICAL  PARAMETERS  ARE  ) 

S?S  jbW'uySfSBaVOSJ'Si&iS&'H8 

211  FORMAT! J9H0THE  INTERVAL  OF  THE  WRITTEN  VALUES  IS  .15) 

END 

SUBROUTINE  GTPSI 

THIS  SUBROUTINE  RETRIEVES  THE  STREAM  FUNCTION  ARRAY  AND  ASSOCIATED 
CONSTANTS  FROM  DISK 

COMMON  PSI ( 121 ,201 1 . UR A T . RC , NPRC , NPZB, R A , NPR , INTVL 

REA0(S)PS1, I TEHS.URAT.UA, RB.NPRB.NPZB.UB.RA.NPZ.NPR.RC.NPRC.NPN 

WR 1 TE (b»  200 ) 

WRITE (b,201 )RA,HB.NPZ, NPH, NPZB, NPRB, UR AT. UA.UB. ITERS, RC, NPRC. NPN 
WR*  TE (b» 202) 

DO  20  Jal.NPZ, INTVL 

20  WR  I  TE tb, 20 3 ) TPSI ( I. Jl. I  a  1. NPR. INTVL) 

200  FORMAT  (20HO  DISK  STORAGE  CHECK)  _  „  „  , 

201  FORMAT !65H0STHEAM  FUNCTION  FOR  FLOW  IN  TWO  CONCENTRIC  PIPES  WITH  P 


FLOW  IN  TWO  CONCENTRIC  PIPES  WITH  P 


IARAMETERS/  10HORA  a  ,F7.2/  10H  HB  a  ,F7,2/  10H  NPZ 

215  /  10H  NPH  a  .15/  10H  NPZB  ■  ,15/  10H  NPRB  a  ,15  / 


3 1  OH  URAT  a  .F9.0/10H  UA  a  .F9.0/10H  UB  a  ,F9.9 

010H  ITERS  a  ,I5/10H  RC  *  .F7.2/10H  NPRC  a  ,I5/10H 

202  FORMAT !05H1 (HE  CURRENT  VALUES  OF  PSI  STORED  ON  DISK  ARE) 

203  FQRMAT(»0'.1SF7„0) 

RETURN 

END 

FUNCTION  CORE (RE ) 

THIS  FUNCTION  COMPUTES  THE  PRODUCT  OF  DRAG  COEFFICIENT 
AND  REYNOLDS  NUMBER  FOR  A  SPHERE  AS  A  FUNCTION  OF 
REYNOLDS  NUMBER 

CONSTANT  COEFFICIENTS 

Alt  1/29. 

A2a-2. 13bJ*l ,E«09 
A3a2. 0159*1 ,£*0b 
A9t-6. 9105*1. £-09 
BOa-1. 2953b 
Bla9.6b*l  E-01 
B2a-9.bb77*l,E-02 
B3«l. 1235*1.6-03 

CHOOSE  THE  APPROPRIATE  POLYNOMIAL 
iF (RE-9. 0)2, 7, 7 
INITIAL  ESTIMATE 

2  IF (RE-0. 00001)3, 9, 9 

3  CORE  a  29.0 


a  ,F9.9  / 
a  , 15/ 1  OH  NPN 


GO  TO  30 
9  X  *29  .  *RE 


BEGIN  NEWTON  METHOD  ITERATION 


k 


0  6  I TERal t 20 


l 

xSxioELx 

CHECK  FOR  CONVERGENCE 


X«A1 #XtA2*X*A2jA J*X**I*A4*X**A»RE 
£XMjj2^*A2  •X*J.*AS*X**2  +  <».*A4*X**J 


tm 


i.X/X>-EP|)5#*jfc 


INITIAL  ESTIMATE 

7  l&y&sBnsjHc’si1" 

BEGIN  NEKTON  METHOD  ITERATION 


?S«SoA^{*5tii*X**a*BjAKAAj  -  ALOO (RE) AELOG 


mi &p* 

CHECK  FOR  CONVERGENCE 

d2lX/X)-EPS)22,22.2<I 
.  **X/RE 

J 

(if 202) 

FORMATS  FOR  OUTPUT  STATEMENTS 
202  FORMAT ( 1  AMO  NO  CONVERGENCE) 


Subroutine  ssm22(n# v#f# x,m,irunb.m) 

waa«C!W!:«^.f,M, 
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PASS  I 

mtv 
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